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Abstract 

We consider off-diagonal Jacobi matrices J with (faster-than-exponential) sparse pertur- 
bations. We prove (Theorem 13. ip that the Fourier transform |/| 2 dp(t) of the spectral measure 
p of J, whose sparse perturbations are at least separated by a distance exp (cj (In j) 2 ) /S J , for 
some c > 1/2, < 5 < 1 and for a dense subset of Cq°(—2, 2)-functions /, decays as i~ 1//2 f2(t), 
uniformly in the spectrum [—2, 2], Q(t) increasing less rapidly than any positive power of t, im- 
proving earlier results obtained by Simon (Commun. Math. Phys. 179, 713-722 (1996)) and 
by Krutikov-Remling (Commun. Math. Phys. 223, 509-532 (2001)) for Schrodinger opera- 
tors with sparse potential that increases as fast as exponential-of-exponential. Applications 
to the spectrum of the Kronecker sum of two (or more) copies of the model are given. 



1 Introduction 



The present paper deals with the Kronecker sum (J denotes the identity matrix) 

K = J®I + I®J (1.1) 
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which are sparse perturbations of the free Jacobi matrix J in the sense that the sequence (p n ) n >o 
differs from the unit on a lacunary subset of natural numbers A = {a^ }j>i, i. e., 



Pr, 



p if n G A , 
1 otherwise , 



;i.3) 



for some p G (0, 1). We assume that A is possibly a random set, df- = cij + Uj, with 



and cij satisfying the "sparseness" condition lim^oo ctj/cij-i = (3 > 1 ((3 = oo is included). A 
concrete example is given by 



dj — Ojj — l 



J 



1,2, 



:i-4) 



with ao = 0, > 3 and ooj, j > 1, independent random variables defined on a probability space 
(fl,B, /x), uniformly distributed over the sets Aj = {—j, . . . These variables introduce uncer- 
tainty in the position of the points where the p n differ from the unit and their support increases 
linearly with the index j (see Remark 1.4 of |CMWlj for less restrictive examples). A disordered 
potential of this type was introduced by Zlatos 0. We shall however consider the deterministic 
case with (II. 4p replaced by a sequence C#j)j>i that increases faster-than-exponential. 

The Jacobi matrices (11. 2p when applied to a vector u = (m„)„>o G ^(^+) can be written as a 
difference equation 

(Ju) n =p n U n+ i +Pn-\U n -l , 

for n > with u_i = 0. We denote by the Jacobi matrix J which satisfies 0-boundary condition 
at —1: 

U-i cos0 — Mo sin = . (1.5) 

is a (noncompact) perturbation of the free Jacobi matrix Jq , where p n = 1 for all n > — 1: 
= Jq + V, the "potential" V composed by infinitely many random barriers whose separations 

increase, at least, exponentially fast. The Jacobi matrix J corresponds to the matrix J° satisfying 

O-Dirichlet boundary condition at —1. 

There have been a few results on models that exhibit spectral transition, supporting spectra of 
different types in complementary set of parameters. The Anderson model in a Bethe lattice [K] is 
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an example. Our intention is to provide another instance of models whose spectrum contains pure 
point and absolutely continuous nonempty components. For the model (II. ip with random sparse 
potential satisfying (jl,4p . the essential spectrum a ess (K) can be decomposed into continuous a c {K) 
at its center and dense pure point a pp (K) near to the edges. Whether a c (K) has an absolutely 
continuous component still remains unknown. 

A method of study the spectrum of sparse Jacobi matrices J, exploiting the uniform distribution 
of the Priifer angles with fixed energy, has been introduced in [MWGA] . With this method, the 
Hausdorff dimension of the spectral measure of sparse block- Jacobi matrices J ® II + I ® Jl can 
be determined with any degree of precision, provided (3 is large enough and Jl homogeneous (see 
[CMWJ). A sharp spectral transition from singular continuous to pure point spectrum, announced 
by Zlatos for a model J* = Jq + V with diagonal potential V (see Theorem 6.3 of [Z]), has been 
proved in [CMW lJ for Jacobi matrices J with random sparseness using this method. Here we 
address the Fourier-Stieltjes transform dp of the spectral measure p of J and investigate the least 
faster-than-exponential sparse condition required for pointwise decay of dp(t). While we consider 
off-diagonal case for reasons explained in |MWGA] . there will be no difficulties of principle in 
applying our methods to the diagonal case considered in [Z]. 

The layout of the present paper is as follows. Some notions and motivations are presented 
in Section [2] and precise statements, Theorem 13.11 and Corollary 13. 2[ are formulated in Section 
Section [3] also contains the mathematical tools used in Section H] to prove Theorem 13.11 None 
of those tools depends on whether the Priifer angles are uniformly distributed. Lemma 13.31 in 
Subsection 13.21 uses the idea of the proof of Theorem 10.12 in Chap. XII of Zy and the Gevrey 
type estimates of Subsection 13.31 permit integration by parts to be applied an unlimited number of 
times. 



2 Preliminaries 

According to [MWGA] (see e.g. Theorem 4.4), if the angles of Priifer are uniformly distributed, 
A G [—2, 2] belongs to the essential support of the singular continuous spectrum of J provided 

r 

r 1 (2 - 1) 

where r = r(p, A) = 1 + $(p) /(4 — A 2 ) and $(p) = (1 — p 2 ) 2 /p 2 is monotone decreasing function of 
p E (0, 1). The local Hausdorff dimension of p in this case reads 

/ lnr \ , 
a H = max 1 - -— , (2.2) 



ln/3' 

(see Theorem 3.11 of [CMW] ). Note that an = chh(p,P, A) varies from to 1 as p varies from p* 
to 1, where p* is defined by r(p*, A) = (3; as A varies from to ±2, a# varies from a#(p, j3, 0) < 1 
to 0, attained at A* defined by r(p, A*) = (3. 
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Let / : [—2, 2] — > C be a smooth function with supp/ E5 A compact and sufficiently localized 



so that ip — f(J)$o is m h(%<+)- The probability |(0,exp {itJ) 
in its initial state ip can be estimated observing that 



of finding at time t the system 



(■0, exp (it J) ip) 



\f(\)\ 2 e ux dp(\) = \f\'dp(t) 



(2.3) 



where p(I) = \\E(I)S \\ 2 is the spectral measure of the state S Q localized at 0, and E(-) is the 
spectral resolution of J. We shall also denote by p^(I) = \\E(I)ip\\ 2 the spectral measure of an 
state ip = f(J)S and the Fourier-Stieltjes transform of p^ shall also be written as dp^(t). 

Our purpose is primary to prove that (12.31) decays to zero as \t\ goes to infinity. Let «f be the 
supremum of a > such that 

\jfdp(t)\ < C (1 + \t\y a/2 , (2.4) 

with C < oo, holds for any t G K, uniformly in a dense (in L 2 (— 2,2)) set of the /'s, sufficiently 
localized around a point A of the spectrum. It follows by a theorem of Frostman (see e.g. [M]) that 
the so called Fourier dimension of satisfies 

(Fourier dimension and Hausdorff dimension do not agree in general). Measures for which ap > 
and o:f = O-u are, respectively, known as Raychmann measures (see e.g. |Ly| ) and Salem measures, 
after Salem's work [S] on continuous distribution functions which are constant in each interval 
contiguous to a perfect set of Lebesgue measure zero. 

Disorder plays a crucial role in most examples of sets and measures in which the Hausdorff 
dimension and the Fourier dimension are equal (see e.g. [HllSa] and Chap. 17 of [K]). The equality 
an = «j? of dimensions defined by (12. 2 p and (12 .4p is thus expected to be attained for models 
with random sparseness condition (jl.4p . However, there are also singular continuous measures 
constructed by a deterministic method (see e.g. |Kot IF] ) and for the deterministic sparse model 
J (or its diagonal version considered in references [Z | IKR] ) with faster-than-exponential (/3 = oo) 
sparsity, it can actually prove that this property is satisfied with an = otp = 1. 

Pointwise decay 



\f\ 2 dp(t) 



<C\t\ 



-1/2 



lnltl 



(2.5) 



(with \t\ > 2, because of the behavior of In \t\ for \t\ < 1) has been obtained earlier by Simon 
[S] for continuous Schrodinger operators with generic and sufficiently sparse potentials and, by 
using a different method, Krutikov-Remling |KR] have found, for a model similar to the one 



considered here, a resonant set 1Z C K. in which 



\f\ 2 d P (t) 



< C (1 + |tir 1/2+£ holds if t e ft and 



\f\ 2 dp(t) 



< C (1 + \t\) 171 for an arbitrary large but finite m, otherwise. The work of |KR] was 



motivated by the "little control" of reference [S] on the rate with which the barrier separations 
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a, — dj-i have to increase (Simon's method requires aj ~ exp(exp(cj 3 / 2 )) and Krutikov-Remling 
need e.g. aj ~ exp(exp(cj)), c > 0, to satisfy their condition aj < Ca^i for some C > and 
fj, > 0)Jj In the present work, the sparseness condition is less restrictive than those stated in 
references [S] and |KR] . Thanks to the Gevrey type estimates developed in Subsection I3.3[ the 
technique of integration by parts, employed in reference |KR] . has been exploited to its limit and 
( 12.51) (with in the place of In Q(t) increasing less rapidly than any positive power of t) has 

been established for aj — aj-i increasing slightly-faster-than factorial: exp (cj'(lnj) 2 ) /S^ for some 

c > 1/2 and 5 < 1 (see H3.4[) for sparseness improvement: if aj — Oj_i = exp ^— j In /(P for some 

£ > 0, then (12.51) holds with In \t\ replaced by |t| e ). 

Our motivation in considering the Kronecker sum ( II .ip is the same that led Simon to investi- 
gate the pointwise decay (12. 5ft : it comes from the observation that the convolution of two singular 

continuous measures, dp^ * dp^p, may be absolutely continuous. Note that dp^ * dp^(t) = dp^ (t) 
is square integrable and this, according to the well-known folklore result, implies that the cor- 
responding measure is absolutely continuous with respect to Lebesgue measure. For model with 
faster-than-exponential sparseness the spectral measure is singular with respect to Lebesgue and 
has Hausdorff dimension 1, uniformly over the essential spectrum [—2,2] (see Theorem 1.4 of [Z], 
which also applies to the off-diagonal model). The spectrum of K in this case is purely absolutely 
continuous (see Corollary I3.2p . We are, however, interested in spectral transitions and to achieve 
the decay ( 12 .4p depending on the local Hausdorff dimension an, a more sophisticate method that 
exploits the randomness of {a^ }j>i is required. Investigation in this direction will be carried in a 
separate paper. 

As in |KR] , our starting point is the representation of the spectral measure as a weak-star-limit 
of absolutely continuous measures (see also [P]): 



_I/ 2 f(A) 

n^co n J_ 2 \y N (X) - w(X)y N+l (X) | z 



f{X)dp{\) = jim - / f(A)- — — — -d\ (2.6) 



for every continuous function f : [—2,2] — > C. Here, y n = y n (z) denotes the solution of the 
eigenvalue equation 

PnVn+l +Pn-lVn-l = z Un , 

satisfying the initial conditions y~i(z) = and yo(z) = 1, and w(z) = z/2 + i^Jl — z 2 /4 is a 
Herglotz function (maps the upper half-plane H into itself). We shall apply this formula with 
f(A) = |/(A)|V\ 

It turns out that y n {X) is a subordinate solution of Ju = Xu (see details in K'.MW. CMWlj ). 
Let T(N, A) denote the 2x2 transfer matrix associated with the eigenvalue equation Ju = Xu with 



1 The above mentioned decay (1 + |i|) 1 / 2+e holds in ref. [KRj for every t £ 1Z if aj-i < Ca^ 2 is satisfied (i.e. 
for fi = 1/2) and a 3 ~ exp (cexp j) with c > ln2. In ref. [S], this decay holds provided aj ~ exp (C e j 3 / 2 ) for some 
a > 0. 
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A = 2 cosf for some ip G (0, 7r) and define <f>j, for the subsequence (iVj)i>i with Nj = aj + 1, by the 
equation 

|T(iV J ,A)| 2 v^=t-S j 

where \T(Nj,X)\ 2 = T*(Nj, X)T(Nj, A) is a real symmetric unimodular matrix, tj = \\T(Nj, A)|| is 

the spectral norm of T(Nj,X) and = f cos< ^ J Under the assumption that the Prufer angles 

\sin0/ 

are uniformly distributed (satisfied if J is random; see |CMW1] ). t 2 = 0(r J ) and we have (see 
Proposition 3.9 of |CM W] ) for an improved version of Lemma 2.1 of [Z]) 

\y Nj (\) - w(X)y Nj+1 (\)\ 2 = \UT(N j: X)v r \ 2 = O(r^) (2.7) 

for A in the essential support of p. Because the hypothesis of Theorem 8.1 of Last-Simon |LS] is 
verified, the limit 0* of the sequence {4>j)j >1 exists and by the Gilbert-Pearson theory 0* = for 
a. e. ip G [0, 7r] in the essential support of p (see |GPj ). establishing the decay property ( 12. 7\i . 

We also have (with A = 2 cos (p and U = I ^ S ^ ) ) 

11- cos <f J 

\UT{N 0l \)w \ 2 = R) , (2.8) 

where Rj = Rj(f), j = 1,2,..., are the radius of Prufer associated with J. Hence, our next 
ingredient is related with the following identity 

\_ I , g (2 . 9) 
for some conveniently chosen j* < j, together with 



-K^ +1 a + ocos20fc+i + csinzfc'fc+i 

where = is the fc-th angle of Prufer; a, 6 and c are functions of p and <£> defined in [MWGA] 

and # = - / H(ip, 9)d9 = 0. Plugging (E9]) and f l2~TU|) into (l2lj) . yields 

io 



|/| 2 dp(0 = - r|/(2cos^)| 2 ^e 2 ^ os ^/ r 

7T ./n ^|*-1 



oo . 2 

+ E - / l/(2cos^)| 2 ^iJ(^,4 + i)e 2 ^ os ^. (2.11) 



Expanding H(ip,9) in Fourier series (see [KR], for details) 



H{ip,9) = {MvTe 2 ™ 6 + A(<p) n e- 2ine ) 

71=1 
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with r given by (12 .ip and tan 5 = c/6p each integral involved in the sum (12. lip is of the form 

/oo 
e ith ^d (G o A) (<f) (2.12) 
-oo 

for some function dG/dX with support 2 cos[y?_, ip + ] C (—2,2), and h of the type: 

hk(t, n; ip) = 2(cos if H — S k (if)) , n G Z 

for some k > j*. 

The main contribution to ( 12.12ft comes from the stationary phase 

2n 

/ i / fc (^) = -2sin^ + T ^(^) = 

and the (^'s that satisfy the equation will be called resonant or critical values. Although the usual 
method of stationary phase does not apply, its ideas can, nevertheless, be easily traced. To each 
integral for which there are no resonant values inside the interval [<£>_, ip + ], the standard stationary 
phase method applies integration by parts once. In order to apply integration by parts m times, 
with m a large fixed number, obtaining therefore a pointwise decay t~ m , Krutikov-Remling in |KR] 
have exploited the fact that dG/dX is analytic and the sequence (aj)j yi were super-exponentially 
sparse. 

We have in the present paper extended Krutikov-Remling's method in two directions. Firstly, 
for pointwise decay t~ l fixed to each of those no resonant integrals, we use infinitely many integra- 
tion by parts to obtain the least possible sparseness condition. Secondly, concerning the integrals 
in which there are resonant values inside the interval [if-, ip + ], we apply Lemma [3.31 Following the 
ideas in the proof of Theorem 10.12 of Zygmund's book on trigonometric series, we use Plancherel 
identity in order to obtain the main contribution of the stationary phase, up to a logarithmic 
correction, provided the integral with the phase e tth ^ replaced by e %tx ^ decays as t~ x . 



3 Faster— than— exponential Sparse Models 
3.1 Statement of Results 

We devote this section to show that \f\ 2 dp(t) decays as t~ 1//2 f2(t), uniformly with respect to any 
Cg° function / with support contained into the essential spectrum [—2, 2] of J, provided p is the 
spectral measure of J given by (ll.2p and ( 11. 3p with the sparseness increment 

%-%-i = /%, J = 1,2,... (3.1) 

2 Note that |A(<^)| < a < 1 uniformly in each compact set K of (0, 7r) and the series is uniformly and absolutely 
convergent. 
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an increasing faster-than-exponential sequence: 

Jim pj-^pj = 0. (3.2) 



From now on, J will always be such that f)3.2p is satisfied. Our goal is to find the least increasing 
sequence that leads to this result. Without loss of generality, we assume that 



< 5 (3.3) 



holds uniformly in j for some < 5 < 1 satisfying^ 

p 2 1 

sup sup < - (3.4) 

^esupp/oA e a + b cos 29 + c sin 29 b 

Note that 

exponentially fast, as j goes to infinity, uniformly in supp/ o A. We now state our result. 

Theorem 3.1 Let p the spectral measure of J associated with the state 5q localized at and let f 
be a smooth function with compact support inside (—2,2) and such that ^ suppf . Suppose that 
the sequence (aj)j>i satisfies 113. l)) -( f5^] ) with 

Pj = ^" ex P (cj (In j) 2 ) (3.5) 

/or some c > 1/2, as j tends to infinity. Then, there exist a constant C, depending on f and p (the 
intensity parameter of J ), such that 



\f\ 2 dp{t) 



< C\t\- 1/2 Q.{\t\) (3.6) 



holds for \t\ > 2, where Q(t) increases less rapidly than any positive power oft. 
Moreover, if the sparseness increments (Pj) ■ >1 is chosen as 

P j = ^ J exv(e- 1 jlnj) (3.7) 

for some e > small enough, then the conclusion ^3. (A) holds with the upper bound replaced by 
C(l + \t\)- l ' 2+£ . 

The proof of Theorem 13.11 uses a classical result on the decay of Fourier-Stieltjes coefficients 
c n (dG) of a monotone increasing singular continuous function G originated from a Riesz product 



3 By equation (4.15) of |MWGA| 6 > mm v -< v < v+ (a/p - \/{a/p) 2 - l) > if < p < 1. 
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(see Theorem 10.12 in Chap. XII of |Zy] ) . The Lemma stated below reduces the resonant estimate 
to a non-resonant one, which will be studied in Subsection |H Our proof resembles, in this sense, 
the proof of Simon (compare Lemma 13.31 below with Lemma 4.2 of [S]) with the non-resonant 
estimate given by Krutikov-Remling's method of integration by parts. 

A well-known folklore result (see, e.g., [C], Exercise 11, Section 6.2) states that if the Fourier- 
Stieltjes transform of a finite Borel measure is square-integrable, the corresponding measure is 
absolutely continuous with respect to Lebesgue measure. We have, thus, as a direct corollary of 
Theorem I3.lt 

Corollary 3.2 The spectrum of K = J ® I + I ® J , with J defined as in Theorem \3.1\ is purely 
absolutely continuous. 

An alternative proof of a theorem that includes the above mentioned folklore result (which has 
been used as early as 1958 to provide nontrivial examples of the statement that the convolution 
of two singular measures may be absolutely continuous; see [KS] ) may be found in |CMW2] (see 
also [SJ, Corollary 3.2). Note that the spectral measure of K associated with the tensor product of 
two vectors in ^(^+) is given by the convolution of the two measures associated to each of these 
vectors (see [SJ for details). 

Of course the corollary extends to the Kronecker sum of any number of copies higher than two, 
and thus the nature of the spectrum changes dramatically in dimension two or higher. This will 
be exploited further in |CMW2j . 



3.2 Basic Lemma 



The goal of this subsection is to prove the following 

Lemma 3.3 Suppose G : R — > M. is a monotone increasing continuous function, with dG sup- 
ported in some closed interval [a,b] C (—2,0) U (0,2), whose Fourier-Stieltjes transform satisfies 

e itx dG(X) < -— J -r- (3.8) 

oo i ~t~ m 

for some constant C < oo and every t G R. Let 7 : R — > C be defined by 

/oo 
e itx ^ x) dG(\) 
-00 

where 

tx(t, A) = tX + -cos -1 - (3.9) 

7T 2 

is a mapping from [—2t, 2t] into \—2t + k, 2t}. If k— K(t) = O (\t\), then 

l7(*)l<-pl*|t| 
holds for some B < 00 and every t el with \t\ > 2. 
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Proof. Denoting by x the characteristic function of the interval [a, b\: x(A) = lifa<A<& and 
= otherwise, by the Plancherel theorem 



OG 



y(t)= e ltx{t ^x{\)dG{\)= A(t,r)dG(r)dr (3.10) 

J— oo J — oo 

where 

A(t,r) = — / e i(te( *' A)+TA) rfA . (3.11) 
2vr y a 



We apply van der Corput estimates (see e.g. Lemma 4.3, Chap. V of |Zy| ) in order to obtain the 
asymptotic behavior of A(t, r) for large t and r. The integral (13. lip is of the form 



rb 

27rA(t,r) = / e 2wifW d\ 

J a 

2tt/(A) = (t + ^A + ^cos- 1 ^ 

7T 2 

where the second derivative /" of / is strictly negative (under the hypothesis ^ [a, b], cos -1 A/2 is 
strictly concave for every A G [a, b}) and proportional to k = O (t), uniformly in r: |/"(A)| > p \k\ 
where 27r 2 p = sup Ae r a6 i A/(4 — A 2 ) 3 / 2 . By van der Corput's lemma 

|A(( ' T)|£ ^ £ 7!TW (3 ' 12) 



holds for a constant K, independent of t and r. On the other hand, for |r| large compared with 
\t\, let us say |r| > A \t\ for A = 1 + |/c/t| sup Ae u b i 1 / (jryA — A 2 ) , the derivative /' of / is of order 
t and Van der Corput's lemma gives 

|A(t,r)|<^ (3.13) 
|t| 

for some constant K', independent of t and r. 

Estimates ( 137121) and ( jSHSj) together with (ETBjl . yield 



r|<A|t| 



| 7 (*)| < / |A(t,r)| dG(r) rfr+ / |A(t,r)| <ZG( 



M>A|( 



r 



# y c J r k' c , 

VI + \t\ J\r\<A\t\ 1 + \r\ J\r\>A\t\ \r\ 1 + \r\ 
< B'\t\~ 1/2 ln\t\ + B"\t\~ 1 (3.14) 
for \t\ > 2 and some finite constants B' and B", concluding the proof of the lemma. 

□ 
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Remark 3.4 If the pointwise behavior of A3. 8]) is replaced by C/(l + |t|) 1_e for some e > 0, then 
the logarithmic correction In 1 1 1 in the conclusion of Lemma has to be replaced by \t\ £ (see the 
last inequality of We shall also apply Lemma \3.3\ with n*K = 0(t) where n* = n*(t) is 

a nonvanishing integer valued, piecewise monotone function, increasing slower than any positive 
power of t. In this case, logarithmic correction will be replaced by yVn*(t) (see Section [^} for 
details). 



Remark 3.5 The function x = x(X) in Theorem 10.12 of Chap. XII of \Zy^ , whose proof has 
suggested us Lemma \3.3\ is a one-to-one mapping of [— it, vr] onto itself (it does not depend on 
t, as in our case; see A3. 9]) ) and j(t) = dF(t), thereby, where F(x) is an increasing function 
with F o x(X) = \/2tvG(\). These facts are not necessary for the conclusion of Lemma \3. 3\ The 
hypothesis on dG(t) in that theorem is, in addition, stronger than ours. Because dG(t) decays only 
"on the average" it is necessary an improved estimate K' /r 2 instead of a simpler one A3. 13]) . 



3.3 Gevrey Type Estimates 

To prove Theorem 13.11 we need to improve results of [ KR] in order to bring (3j down from the 
asymptote exp (exp cj) to (j!) 1 ^ for some e > 0. The following proposition gives a Gevrey type 
estimate for the derivatives of the Priifer angles. 



Proposition 3.6 Let 9^ = 9k((p), k = 0, 1, . . be the sequence of Priifer angles starting from 9q, 
consistent with the initial condition ( 2.7 ). Then, for every m > 1, 

\9' m (cp)\ <C x 8 n p m (3.15) 

and for every m > 1 and n > 1 

-Ad^(<P)\<C n V n ^-i (3.16) 
hold uniformly in compact subsets of (0, ir) with C n = K/n 2 , K < 3/(2tt 2 ) = 0.151981... and 

1 + A 



V 



5K 



where A < 1 (A = 0(5) ) is a constant satisfying 9' m (<p) < (1 + A) (3 m , which is computable from 
Proposition 5.2 of IMWGAf . and 5 is suitably small. 



Remark 3.7 Proposition ^. 6\ replaces the unspecified constant Cj appearing in Lemma 3.1 of [KR] 
by Krfj\/j 2 . Detailed information on the growth of 9m \<p) in both j and m are an essential 
ingredient of our method. 



To prove Proposition 13.61 we need the following 
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Lemma 3.8 Let C * D denote the convolution product in M. z+ : 

n 

(C*-D) n = J2C t D n . l (3.17) 



i=0 
;2 



for n > . If C has components given by C{ = K/i , i > 1, with Cq = K , then 

(C*C*---*Q < C 

k— factors 

holds for every k > 1 provided K < 1/ (2 + 27r 2 /3). // Co = 0, iaen same result holds with 
K < 3/(2tt 2 ). 

Remark 3.9 The n-th component of the convolution of C with itself k -times satisfies 

C n ---C lk <C n (3.18) 



h,...,i k >0 
U H Hfc=n 



6?/ Lemma \3. #1 



Remark 3.10 VFe aave stated Lemma \3. 8\ for sequences C = (Cq, Cx, . . .) G M z+ u^a tae O-ta 
component Cq = and Cq = K since both cases will be considered in this subsection (see Proposition 
\3.13\ below for the case Cq ^ 0). Lemma \3. 8\ plays a key role in every estimate involving higher 
order chain rule. 

Proof of lemma. By definition, 
Writing 

n 2 _ n 2 (1 1 \ 

i\n-i) 2 ~ ( n - l ) 2 + l 2 {~ 2+ Jn~^f) ' (3 ' 20) 

the pre-factor in the r.h.s. of (I3,20p can be bounded using < (a - b) = 2 (a 2 + 6 2 ) - (a + 6) 2 , 
which holds for any real numbers a and 6, with a = n — i and b = i. We have 



< (n - 2if = 2 ((n - i) 2 + i 2 ) - 

or, equivalently, 



u 2 



n 2 



fn — i) 2 + i 2 



< 2 . (3.21) 
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Plugging fl3T20|) into fl3TT9|) together with ( KTflf . gives 



^ 2 



< 2iT (^1 + —j < 1 (3.22) 

provided K < 1/ (2 + 27r 2 /3). The case Cq = 0, the terms with i = and n do not contribute to 
the sum and the inequality ( 13 .221) holds provided K < 3/(27r 2 ). Once we have C * C < C, Lemma 
13.81 is proved by induction. 

□ 

Proof of Proposition \3.b\ The proof uses the recursive relation 

m {f) 

— go 6 m -i(ip) + (3 m ip (3.23) 

where 

g = g ((p i 6) = tan -1 ((tan 9 + cot (f) /p 2 - cot if) (3.24) 
together with the Scott's formula for higher order chain rule (see e.g. |FLy| ) 

n 

{9°ft ] = Y.9 [k] °f E f ln] ---f [lk] (3-25) 

k=l ,ife>l 
ilH Mfc=rt 

where, from here on, stands for h^/nl, the n-th derivative of h divided by n\. 
Upper and lower bounds for the first derivative has been provided in |M WGA] : 

(1 - A) m < Q'Jif) < (1 + A) f3 m (3.26) 

with A < 1 a constant. Now, choosing r] = — ^— , (13.261) establishes ( 13 . 1 5[) for every m > 1. 

oK 

Since • >1 is a fast increasing sequence we apply the Scott's formula to g in ( I3.24[) as it were 
a function of a single variable 9. This really gives the main contribution to the derivatives, g as a 
function of z = e t6 , continued to the complex plane, is analytic outside a disc of radius strictly less 
than 1 — e/£ < 1. The derivatives of q(e l6 ) = g{9) may be estimate by Cauchy formula: 

\g [k \0)\ < Ci e (3.27) 

holds for k > 1 with c± as small as one wishes, by increasing £ accordingly ( (13.271) can be bounded, 
e. g., by e (ci^/e) k = e£ k , for any e > 0). Replacing / by Q m -\ in (13.251) . gives 

n 

0t\v) = 5> W ° ^-iM E • • • ( 3 - 28 ) 

fe=l «i,...,ifc>l 
ilH hifc=n 
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We prove (I3.16P by induction in n. Consider the case n = 2, for any m > 1. By equation (I3.23p . 
together with (13.271) and (13.151) . we have 

C = g"o9 m ^-{9' m _ 1 f + g'oO m ^-9l_ 1 

The iteration of this relation together with ( I3.27P and (13. 3p . yields 

i=i 



m-l q2 

m-j 



Pm-l 



1 - c^S 2 



< 2C 2V 2 (3 2 m _ 1 
provided S is chosen so small that Ci^5 2 < 1 and 

are both satisfied, establishing (I3.16P for n = 2. 
Now, suppose 

holds for m > 1 and j = 2, . . . , n — 1 and we shall establish the inequality for n. By this assumption 
together with (I3.3p . we have 

€ ] -i(v?) < (W >!„ , = erf ( ^ < r ;( ^, , . (3.29) 



Plugging ( 13TT5|) . ( l3~2Tj) and ( 13T29|) into ( BHD , together with (l3TT8|) . yields 

n 

fe=2 ii,...,ife>l 

ilH h«fc=n 

Here, we have separated the term with k = 1 which applies n derivatives on Q m -\- Note that, for 
all the other terms with k > 2, we have ii, . . . , > 1 and the derivatives applied on the 9 m -\ are 
of order strictly smaller than n. The iteration of this relation, gives 

C rn-l nn 

Q{n]f,„\ „ S r i (xc^\riQn ( n fV'-l ' m-j 



? 1 =1 Pm-l 
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provided 

c iT^T -i ^TT- < 1 ( 3 - 3 °) 

holds for every n > 2. We pick 5 satisfying both ( 13. 3 p and (I3.30p for n > 2, concluding the proof 
of Proposition 13.61 

□ 

Remark 3.11 The well known formula for higher derivative of composite functions, Fad di Bruno's 
formula, cannot be used recursively since the constant rj in equation A3. 16]) deteriorates each time 
it is applied (see eq. (6.10) in Subsection 6.2 of IBMtf ). The proof of A3. 16]) . by induction, using 
Scott's formula was based on yet unpublished manuscript "O(N) Hierarchical Model Approached by 
the Implicit Function Theorem" by W. R. P. Conti and D. H. U. Marchetti. 

To obtain the t _1 decay from the summation in (12. lip , it is necessary to apply an arbitrarily 
large number of the integration by parts for integrals of the type (12.121) : 

ithj.„ • / / f \ „ith 







f e im dip = i —Jo e^dip 



where supp/o = C (0, 7r). The following propositions gather tools to implement the 

estimate. 

Proposition 3.12 Let fo(ip) and g(ip) be, respectively, C°° complex and real-valued functions on 

[0, 7i") and let L = — g(ip) be an operator defined in this space. If 
dip 

, 1 rn , 1 d d d 
fn = —Lf Q = ——Q—Q---—gfo 

n\ n! dtp dtp dip 

denotes the n-th application of L over fo, divided by factorial of n, n — 0, 1, . . ., then 

/»= £ g^.-.g^f^ . (3.31) 

fcl,...,fe n ,p n >0 

fclH \-k„+p„=n 

Proof. The proof is by induction. For n — 1, we have /i = f//o + Qf' Q . Assuming that (13.311) holds, 

fn+i = — — r (g7n + 
n + 1 

' £ (h + ■ ■ ■ + k n+1 + p n+l )Q^ ■ ■ ■ gi^]f^] 



n + 1 

fci,...,fc„+i,p„+i>0 
fclH hfc„+i+p„+i=n+l 



V" . . . gfcn+l] flPn+l] 



fci,...,fc n+ i,p n+ i>0 
fciH hfc„+i+]3„+i=n+l 
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where we have applied the product rule 

n 

r„ = E E ^•••(^ ] ) / --^ [fei] /o w + E ■ ■ ■ e [kn] (ft ] ) , 

J = 1 ki,...,k n ,p n >0 fci,...,fc„,p„>0 
feiH hfc,i+Pn=" fciH \-k„+p n =n 

used (V^ 1 )' = (fcj + for j = 1, . . . ,n (analogously for (/o Pn] ) = iPn + l)g [Pn+1] ), redefined 

variables: 

k j+ i = kj + 1 , fci+i = fc| for / ^ j and p n+1 = p n 

(the same for j3 n+1 = p n + 1 and = kj for j = 1, , . . . , n) and have added a new variable k 1 . 
Observe that kj+i = kj + 1 > 1 but we can start the sum over kj + i from since is 
identically at kj+i = 0. This completes the proof of the proposition. 

□ 

1 sin 2 
In our application, g = — — — and fo = 1/(2 cos (p) I 2 — . A n (ip) (or its complex conjugate). 

The main contribution for m > j* + 1, where j* = j*(t) is such that 

f3j* < t < f3 3 * +1 , (3.32) 
comes from the derivatives of the Priifer angles and in this case it is thus sufficient to consider 

For m < j*, we have 

^m+i(^) = -2tsin^ + n^ +1 = -2t sin <p (1 + 0(1)) (3.33) 

and the derivatives of higher order 

M M „ fl W _ / 2t(-l) (fc+1)/2 sin^ if fcisodd 

I 2r(— 1) cosy? it /e is even 

satisfies, in view of ( 13.1 6p . 

thl ) +1 ^)<2\t\+nC k r 1 k (3 k m k\ . (3.34) 
It is also sufficient to consider in both cases 

A 



mil 



Ro J- \ a + & cos 2% + c sin 20j 



j=i 
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where F{6) satisfies, by direct computation, 



F ik]<tilK(t\ k (335) 



8 \5 / 

for k > and some positive number (. We also need 

8^x) = t£,k = l,... . (3.36) 

Proposition 3.13 Let 9k = 6k{(p), k = 0, 1, . . be the sequence of Priifer angles and let r], 5 and 
{C n } be the constants that appear in Proposition ^ .b\ Then, there exist positive numbers d and fj, 
which can be expressed in terms of the previous constants, such that (with g = g^ < d/(3 m+ i) 

Q [n] = {soe' m+1 ) [n] <-^-C n r^ (3-37) 

Pm+l 



as well as (f = ff l < Rf5~ m ) 



(m \M 

f ™ = l%\J[ Foei ) -M CniCv)n ^ (3 - 38) 



hold for every non-negative integer n, with ( as in A3. 35\) . 



Proof. These inequalities are established as in Proposition I3.6[ by using the Scott's formula. We 
begin with ( 13.371) . If fj is the smallest constant such that (i — l) 2 ^ 1 /i < fj 1 holds for every i > 1, by 
(13.16!) . (I3.29j) and (13. 25ft with g and / replaced by s and 6' m , we have 



n 

e w = ^/ic^ fa + wK 1 • • • (** + w 

k=l ii,...,«fc>l 

»iH Hk=n 



k=l \"m+l) ii,...,ife>l 
iiH Hk=n 



< dC n f] n ^ 

Pm+l 

where fj = 8fj 2 / (1 — A) and d = Sfj/ (Sfj + A — 1). In the third inequality we have used the lower 
bound (I3.26P for 9' m and (I3.18p . Note Sfj > Sr] > 1, by definition of rj in Proposition 13.61 
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For f)3.38p . we start with the Scott's formula (13.251) with g and / replaced by F and 9j which, 
together with flOSJ, ( KWf and (I3733|> - gives 



ll*fcj 
'i 

fc=l ii,...,ife>l 
ilH Mfc=n 



fe=i v ' 



< lc n (Cr/)™/3; . 

Now we take the n-th derivative of the product. For this, we use the variation of Lemma 13.81 
mentioned in Remark 13.101 



\j=l / rai,...,n m >0 

raiH \-n m =n 



\n m \ 

m) 



^ C n(CvT-^P m (3.39) 
concluding the proof of this proposition. 

□ 



Remark 3.14 An estimate of K3.31 ) with g = l/th' m+1 (ip) for m < j* is analogously given by 

[k] - d. 



Q [k] = (soth' m+i y Ki <^C k r) k p k m 



(3.40) 



with d n = 2dn/c where c = mm vesup pf \2sm<p. Note that, for any fixed m and t satisfying A3. 32) , 
the second term of the l.h.s. of 34\ ) rapidly overcomes t. On the other hand, an estimate for k 
in which t still dominates \3. 34\ ) is, by the Scott's formula A3. 25\) . much better than \3.40\j : 



< 



- V oh' V — h {il+1) 

1=1 ii ,...,€,>! l ' 



h,—,ii>l 
ilH h«;=fc 



A 2 l 1 1 

1=1 h,...,ii>l 
»iH M;=fc 



1 



h 



IV- 



I m+1 



< 



2 fc 1 A 2 l l k 

]t\ Jfe! ^ d+i ' 
1 1 i=i 



(3.41) 
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Let us put all together. Plugging (13.371) and f!3.38j) into (13. 3 1 j) . deduced in Proposition 13. 121 by 
applying n times integration by parts n\f n = L n f to the integrand fo of (I2.12p . we arrive at the 
following estimate: 

Jkl] Jkn] f\Pn] 

' ' ' y JO 



n\ \ f n \ < n\ 



E 



ki,...,k n ,p n >0 
fciH ik n +p n =n 



< -4f n — 

- Rl 5 m 



flu 



~ Rl 5" 



Pm+l 
fl 



nl 



E 



ki,...,k„,p n >0 
fciH \-k„+p n =n 



IV. 



... (3.42) 

L u \Pm+l/ 

where D = d- max (?), (r/). Estimate (13.421) will be used to get an upper bound for all non-resonant 
integrals of fl2~lH . 

Remark 3.15 As C n = K/n 2 , n > 1, (C = K) with K < 1/ (2 + 27r 2 /3) are bounded constants, 
jjj3.4S\ ) makes explicit the dependence on the number n of times that integration by parts is applied 
to integral of type 112. lty) . Explicit dependence of n was not necessary in reference IKRf . since n is 
an arbitrarily large but fixed number. Apart this, {3.^2$ agrees with the estimate used on p. 522 of 

mi- 



4 Proof of Theorem 3.1 



Let t be a fixed number. We assume t positive but the negative value can be dealt similarly. Let / 
be a C°° function with compact support in (0,2) and let If = be smallest closed interval 

that contains supp/o A, X(ip) = 2 cosy?. Since the spectral measure is symmetric, dp(—X) = dp(X), 
we need only to consider / supported in one-half of the essential spectrum. We have excluded 
the origin to avoid that the curvature of cos -1 A/2 in (13.91) vanishes (see observation right before 

dsn). 

For j* = j*(t) defined by equation (13.321) . let n* = n*(t) be given by 

(n* - l)/3j* < t < n*/3 r . (4.1) 

Since f3j*/t > 1/n* there are at most n* points <fi, . . . , (p n * in the support of / o A satisfying 

W' 

-sin^ + -^ = 0. (4.2) 

Observe that 1 < n*(t) < (3j*+\//3j* + 1, by (I3.32f) and (14. lft . For the sparseness increment (3j in 
( 13. 5p . we have 

. i . lnt / „ /In In In 



cln Int \ \ In In t 



19 



and, consequently, the number of resonant values n* is a monotone nondecreasing function of t in 
each interval (ftj*, Let L(t) denote the continuous interpolation of n*(t). It follows from 

these observations that L is a piecewise linear function with inclination 1 / /3j satisfying 

1 < L(t) < Ee cln2lnt = n 2 (t) , j <t</3 j+1 (4.3) 

for some constant E, independent of t. The inequality (14. 3 ft will be used at the end of this section. 
By (12.111) . the Fourier-Stieltjes transform of p can be written as 

oo oo 

\f\ 2 dp(t) = I r _ lfi (t) + Y,(lj,n(t)+ Ij,n(-t)) 

j=j*-l n=l 

where 

1 /*7T • 2 

hn{t) = ~ / l/(2cos^)| 2 S -^A n (<p)e 2it ^ +ne ^d<f . (4.4) 
71 Jo -K-j 

We apply integration by parts to all terms of this sum not satisfying the resonant condition 
(14.21) . Since the support of / is compact, every boundary term vanishes. Integration by parts may 
be repeated Nj times depending on the index j of the sum. Propositions 13.121 and 13.131 together 
with its combined estimate (I3.42p . can be used to get an upper bound for each integral (I4.4p with 
j > j*. This yields 

oooo oo / oo .. \ .. 1 , n \ Nj 

j=j* n=l 3=3* \n=l / 

where a = sup^ e// |^(y?)| < 1- If the sequences (/3j)j and {Nj). are chosen so that 

1 / 8- \ Nj 1 

lAT, ■ 1 / > J 3 \ at i ^ 1 



C M^J (46) 

then the series in 1 14. 5ft converges uniformly in t. By ( 13. 31) and (13.321) . we have 

1 1 Pj. +1 _ _ _ _/^_ < 5j _ r 1 



Pj+l Pj*+1 Pj*+2 Pj+l Pj*+1 

and t//3j* + i < 1. Consequently, 



OO OO f-. oo 1 

EE(M + IW-«)I) * ^E c ».^ 

j=j* n=l u j=j* JT 



< JfLJ I «^f : ,.<£ (4.7) 



holds with C < oo independent of £. 
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Let us now verify that the sparseness condition stated in Theorem 13.11 satisfies (14. 6p . Choosing 
as given by ( 13. 5 j) and Nj = j + 1, by the Stirling formula, 



^fe) % ^vra(^r mp( - 2co+i)in(j+i H i+ KT)) 



and (14.61) holds for any c > 1/2 provided j is large enough. 

Note that, for the sparseness increment 0%L>i given by (13. 7p . 

A', 



(J^) ' = 5 £U+l) exp (-(j + 1) In (j + 1) - e- 1 ^ + 1)) < l^D' 



N, 



and (14. 6 p holds with l//3j +1 replaced by 1/(3^1, provided 1 + 1/e > InD — (1 — e)\n5 and j is large 
enough. Together with Remark 13 A\ the proof may continued exactly as for decaying t~ x . We shall 
consider only the latter case. 

For j < j*, (I3.33P holds and we need replace the estimate (I3.37P by ( I3.40p and t occupies now 
the place of (3 m +i m (I3.42j) . Applying successive integration by parts to If-i,o gives, analogously 



for some constant C Note that, by ( I3.32p . 



and by (jMD 

1 P N * 



D Nk — —± NJ < 1 (4 9) 

Pk+i 

for /c large enough. 

oo 

It remains to estimate the sum Sj*(t) = (lj*-i jn (t) + 7j*_i )Tl (— t)) which contains the most 

n=l 

significant terms responsible for t -1 / 2 decaying behavior. To extract this decay we write 

Ir-iAt) = - r l/(2cos^)| 2 S ^B^_Me 2U ^ +n ^ t U ( p (4.10) 
71 Jo 

(analogously for Z,»„ ln (t)) where, by (13 .23 p . B k is a function of the Priifer angles 9k(f) such that 
\Bk\ = \A\ and 

arg B k = arg A + 9 k+1 (<p) - /3 k +i<P 
= &rgA + g o9 k (ip) 
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with g given by (13.241) . We then apply Lemma [3731 to (I4.10p with 

d(G n o\) fo,) = - |/(2 cos iptf^-B^^dtp 



and 



tx(t, A) = tA + 2n/3j* cos 1 ^ . 



Note that, by (13.321) . k = 2nn(3j* and n* /3j* = 0(t). In order to fulfill all assumptions of Lemma 
it remains to show that dG n {t) decays as (see equation (13.81) ). 

We estimate the Fourier-Stieltjes transform 



dG n (t) 



71 



|/(2 cos p)| 



2 Sin if Dll ; x _2itcosy, 



-fl-»_ 1 



5™_ 1 (^)e 2itcos ^^ 



(4.11) 



as the non-resonant integrals (14.41) with j < j* (see Remark l3.14p . The estimate (I3.37P is replaced by 

sin 2 

(I3.4ip and /o = |/(2cos<^)| 2 — ^ -^-A n (<^) exp (ing o 8j*-i(ip)) includes now an extra exponential 

•\r |U") 

term depending on 9j*-i(cp). 

To deal with this new term we need some more estimates. By the Scott's formula (13.251) together 
with (I3.27p . we have 



(e in9 ) 



[N] 



N u 

Err 

k=l ii,...,i k >\ 
h+-+i k =N 



< -,C V (4.12) 
applied once again together with Proposition 13. 6} 



with C2 = e nCl — 1. The Scott's formula 
( 13729]) and iil2| . yield 



N 



fe=i 



ii,...,i fe >l 
ii+-+ifc=JV 



with C3 = C2^/(£ — 1). Finally, we shall replace f !3.39|) by 



[JV] 



ni,...,^* >0 
TllH hn,'* =JV 
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with ( = max (<5£, (). 

Since the require estimates didn't change significantly, we integrate (14. lip by parts iVj*_i times 
and use (E32D with (EZID in the place of ( 13371) to get, exactly as for (Oil . 

C" 



dG n (t) 



< — d k ~a n 
~ t 



0k <t < f3 k +i 



for some constants C" < oo, d < 1 and a = e ci sup^gj (^(v 9 )! < 1> as c i is arbitrarily small by the 
observation after (I3.27p . Here, we have used the fact that (14. 9 p holds with D replaced by D/d, for 
any d > 1, provided k is large enough. This immediately imply, by a slight modification of Lemma 
EH (see Remark E3D, 

< / |A(*,r)| dG(r) dr + 0(l/t) 

J\r\<A\t\ 

< ^d n Yd k / —dr + 0(l/t) 
^ fo hu T 



< 



2KC" ...^^ Wi 

fe=0 



:n*~a n 



X> fe ln 



0(1/0 



and by (14. ip and the fact that ln/3fc+i//3fc = 0(ln 2 fc), we have 

oo 

|S,.(t)| < E (l^-l,nWI + |^-l,nH)|) 



< 



n=l 
C 



where Q(t) is defined in (14.31) . 

Now, we show that Q(t) increases slower than t £ , for any e > 0. Suppose, by contradiction, that 
lim^oo Q(t)/t £ = k > holds for some e > 0. Then, by L'Hospital, 

,. n(t) ft'(t) c ft(t) ,. lnlnt n 
hm = hm = - hm • hm — = , 

t^roo t £ t^foo st £ 1 £ t^oo t £ t-toc In t 

concluding the proof of Theorem 13.11 



□ 
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